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The standard non-equilibrium theory of noise in ideal bolometers and microcalorimeters fails to 
predict the performance of real devices due to additional effects that become important at low tem- 
perature. In this paper we extend the theory to include the most important of these effects, and find 
that the performance of microcalorimeters operating at 60 mK can be quantitatively predicted. We 
. . . give a simple method for doing the necessary calculations, borrowing the block diagram formalism 

from electronic control theory. 
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Introduction 

<: 

■ A complete non-equilibrium theory for the noise in simple bolometers with ideal resistive thermometers was given 
(N| : by J. C. Mather in 1982 [3 and extended to microcalorimeter performance two years later Q ■ Here we use the terms 

bolometer and calorimeter in the conventional sense, respectively indicating power detectors and integrating energy 
^""^ ' detectors. 

This theory shows that the performance of these devices improves dramatically as the operating temperature is 
' reduced. However, at temperatures below ^200 mK, it becomes increasingly difhcult to construct a bolometer that 
\ behaves according to the ideal assumptions. The resistance of the thermometer becomes dependent on readout 
' power as well as temperature and equilibration times between different parts of the detector become significant. 

■ Thermodynamic fluctuations between internal parts are then an additional noise source. The physical description for 
I most of these effects is straightforward, but combining all of them into a detector model can be algebraically daunting. 
' Theoretical models that describe complex thermal architectures are necessary to understand the behavior of real 

(~| , devices, and some groups have already extended the "ideal" model developed by J. C. Mather in 1982 to include 
Q^' some non-ideal effects in order to explain their experimental results |^ H, 1^ - We developed a general bolometer 
and microcalorimeter model using the block diagram formalism of control theory. The formalism helps with the 
mechanics of the problem, while keeping the physical model reasonably transparent 0; S ■ model we have 

^ I included the thermal decoupling between the electron system and the phonon system in the thermometer, the so called 
hot-electron model, the thermal decoupling between the absorber and the thermometer, and non-ohmic behaviors of 
the thermometer in addition to the hot-electron effect. The hot-electron model assumes that the resistance of the 
thermometer depends on the temperature of the electrons, and there is a thermal resistance between the electrons 
and the crystal lattice through which the bias power must flow, increasing the temperature of the electrons above the 
, temperature of the lattice, and therefore changing the thermometer resistance. This effect is well known in metals at 
low temperatures and has recently been studied in semiconductors in the variable range- hopping regime |lOl IllL Il2| . 
The noise analysis incorporates terms for thermometer Johnson and 1// noise, amplifier noise, load resistor Johnson 
noise, and thermodynamic fluctuations between the electron and phonon systems in the thermometer as well as 
between the absorber, the thermometer, and the heat sink. In the model we also included the effect of thermometer 
non-ohmic behavior, i.e. dependence of the thermometer resistance on the bias signal. This effect is particularly 
important when Transition Edge Sensors (TES) are used and makes the model valuable for predicting the 
performance of this type of detector. 
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FIG. 1: Thermal sketch of a bolometer or microcalorimeter 



I. THE IDEAL MODEL 



To help the reader understand the algebra of our model we decided to start our analysis with an overview of the ideal 
model that has been previously developed. Despite their different applications, bolometers and microcalorimeters are 
very similar detectors and the theory of their operation is largely the same. The considerations of this paper apply 
to both kinds of detectors unless otherwise specified and we will use the generic term "detectors" to refer to both. 

Typically a bolometer or a microcalorimeter is composed of three parts: an absorber that converts the incident 
power or energy into a temperature variation, a sensor that reads out the temperature variation, and a thermal link 
between the detector and a heat sink. The sensor is typically a resistor whose resistance strongly depends on the 
temperature around the working point. An ideal detector can be represented by a discrete absorber of heat capacity 
C in contact with the heat sink through a thermal conductivity G (see Fig. and a thermometer always at the 
temperature of the absorber. The thermometer sensitivity is specified by: 

TdR 

where T is the detector temperature and R is the sensor resistance. The thermal conductivity G is defined as: 

where P is the power dissipated into the detector. The conductivity G can generally be expressed as a power law of 
the detector temperature T, i.e. G = Gq ■ . Notice that numerically Gq is equal to the thermal conductivity at 
1 K, but dimensionally Go is a thermal conductivity divided by a temperature to the /3. 

In equilibrium, with no other input power than the Joule power P used to read out the thermometer resistance, 
the equilibrium temperature of the detector T is determined by integrating Eq. |2] between the heat sink temperature 
Ts and the detector temperature: 

G{T')dT' = P{T). (3) 

Ts 



Assuming the power law expression for G introduced before and integrating it becomes: 

{T^+^ -TI+^) = ^P{T). (4) 



It is important to remember, when calculating the equilibrium temperature, that the power P depends on the value 
of the sensor resistance and, as a consequence, it depends on the temperature T, as explicitly indicated in Eq. ^ To 
calculate the equilibrium temperature it is therefore necessary to solve the system of equations represented by Eq. 01 
the P vs. R curve and the R vs. T curve. In general, the system must be solved numerically. 

Of interest from the point of view of the detector operation, is how the temperature rise AT above the equilibrium 
temperature depends on an external incident power W . The power input to the detector [W + P) is partly stored into 
the heat capacity of the detector and partly flows to the heat sink through the thermal conductivity. The equation 
that determines the generic temperature To of the detector is therefore: 

dT fTo 

+ / G{T')dT' + P{Td), (5) 
at Jt. 
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FIG. 2: Typical detector readout circuit. 



where we explicitly indicated that the bias power can be a function of the temperature Td and where the quantities 
Td, W ^ and P can be a function of time t. We can express the generic detector temperature To as a function of the 
equilibrium temperature T defined in Eq. 0]as To — T + AT. Eq. Othen becomes: 



rT pT+AT 
— + / 

dt Jts jt 



^ d{T + AT) ^ I Q^j^^^j.,^ I G{T')dT' = W + P{T + AT). (6) 



If we stay in the so called small signal limit, i.e. we assume that AT is small compared to T, we can expand the 
second integral to lowest order in AT/T, obtaining: 



^diT + AT) 



T 

^+ / 

dt 



[ G{T')dT' + G{T) ■ AT = W + P{T) + AP (7) 
Jts 



with AP = P{T + AT) — P{T). Subtracting Eq. |21from Eq. [3 and considering that the equilibrium temperature T 
does not change with time, we obtain: 

C^^^^ + G-AT = W + AP, (8) 

where for simplicity we expressed G = G{T). 

In general, the bias power will change with temperature, since R changes, and its expression depends on the bias 
source impedance. A typical bias circuit is illustrated in Fig. [3 where R is the thermometer resistance and Rl is 
a load resistor. The most commonly used bias conditions are near current bias {Rl 3> R) and near voltage bias 
{Rl <C R). More complex bias circuits are also used, and can always be represented by the circuit of Fig. [fusing 
Thevenin equivalence theorems. Differentiating the expression for the Joule power P ~ PR = V"^ / R and using the 
bias circuit of Fig. |21we obtain: 

AP ^ ~-?—J^aAT. (9) 
TRl + R ^ ' 

This term is generally referred to as the electro-thermal feedback term and it often plays an important role in the 
response of a detector. For simplicity in the small signal analytical calculations, we write the electro-thermal feedback 
term as: 

AP^-GstfAT, (10) 

where 

_ P R-Rl 

LrETF = Oi, ii) 

TRl + R ^ ' 



In(a)) 



H((D) 



Out(oi)) 



FIG. 3: Block diagram representation of a system with transfer function H{iL 



SO that Eq. IHI becomes: 



C^^^^ + {G + Getf)- ^T^W (12) 
at 



or, introducing an equivalent thermal conductivity Geff — G + Getf (which we refer to as effective thermal conduc- 
tivity) : 

C^^^ + Geff -AT^W. (13) 

The easiest way to solve this differential equation is using Fourier transforms. The procedure is to use Fourier 
transforms to convert the terms of Ea. ll3l to the frequency domain, solve the equation in the frequency domain where 
it becomes a linear equation, then Fourier invert transform the result to the time domain. The advantage of solving 
Eg. 1131 in the frequency domain comes from the fact that the expression dAT{t)/dt in the frequency domain becomes 
jLuAT{uj), where we used the engineering notation j = V—T- EquationElin the frequency domain then becomes: 

jLuCAT{u) + Geff ■ AT{lu) ^ W{uj). (14) 

whose solution is: 

AT{cu) = -^—4 W{u;) (15) 

Geff l+JUJTeff 

with Teff = C I Geff. 

The detector system behaves as a low pass system, with time constant Teff. For negative electrothermal feedback 
Getf must be positive and the detector time constant is shortened. For positive feedback Getf is negative and the 
detector time constant is lengthened and, in the case of \Getf\ bigger than G, the detector becomes unstable. The 
sign of Getf depends on the sign of a and on the bias condition used (i.e., the ratio R/Rl). In the small-signal 
(linear) limit considered here and in absence of amplifier noise, the signal has no effect on the detector performance. 
However positive feedback reduces the effect of amplifier noise, while negative feedback helps linearize the large signal 
gain, and improves microcalorimeter resolution for large signals at high count rate. Since it can usually be arranged 
that amplifier noise is negligible, these practical considerations normally favor negative feedback. People then use 
current bias {R < Rl) for detectors with negative a and voltage bias {R > Rl) for detectors with positive a. 

In operating a detector, what is really detected is not directly the temperature variation AT, but the resistance 
variation Ai?, which is read out either as a voltage or current variation, that is: 

AV = V- AT (16) 

TRl + R ^ ' 



a R 

AI = -/ — AT. (17) 

We can generically indicate the output signal as X and the relation between the output and the temperature as: 

AX . AT 

— =a-Atr— (18) 

where Atr is a dimensionless parameter that quantifies how much the output signal is sensitive to resistance changes 
and that we call the transducer sensitivity. Numerically Atr is defined as: 

RdX 

Notice that the expression of Atr can be easily derived from Eos. 1181 and 1161 or 1171 for voltage and current readout, 
and is always smaller or equal to unity for passive bias circuit {Rl > 0). 
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FIG. 4: Some general operations with the block diagram algebra. 



The response of a detector is usually quantified by the responsivity Sioj), defined as: 

AX{lj) 



W{uj) 



(20) 



that is, the responsivity characterizes the response of the detector, AX , to an input power W. In the ideal model 
just described we can combine Eas.ll5landll8lto obtain: 



AX{u;) 



1 



1 



X -a- At 



Geff 1 + juJT, 



T 



Wico), 



and the responsivity is then equal to: 



Situ) 



1 



X -a- Af 



Geff 1 + juJTeff T 



(21) 



(22) 



A detector at the working point is also often described by the complex dynamic impedance Z{uj) — dV{uj)/dI{uj). 
The dynamic impedance Z{lu) differs from the detector resistance R ^ V/I due to effect of the electro-thermal 
feedback. When the current changes, the power dissipated into the detector changes too, therefore the temperature 
and the detector resistance change. It is often useful to express the detector performance and characteristics in terms 
of the dynamic impedance since it can be easily measured experimentally. 

The calculation of the analytical expression of the dynamic impedance is simple. Differentiating Ohm's law, 
V — I ■ R, we obtain: 



dV = I ■ dR + R- dl. 



(23) 



Using Eq. |S1 in the frequency domain with W ^ Q and the definition of the thermometer sensitivity a in Eq. ^ we 
obtain: 



dR = -adT - 5^ i dP. 

T GT\+ jujT 



(24) 



with T = C/G. This r is called the "intrinsic" or "thermal" time constant of the detector. Differentiating the 
expression of the Joule power dissipated into the thermometer P = y • /, we obtain: 



dP ■ dl + I ■ dV, 



which, combined with Eas. 1231 and 1241 gives: 



dV = R-dl + I 



Ra 1 



-{I ■ dV + V -dl). 



GT1+ jujT 

Notice that in Eas. 1231 through most of the terms are function of the frequency uj. Solving Eq. [5^we obtain: 



Z(c.) 



dl{uj) "^I-PJ^+JUJT \+ JOJT^ 



(25) 



(26) 



(27) 
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FIG. 5: Block diagram representation of a feedback system. 



where we used the expression 



1 + ^ 

Zo = Z{uj = 0) = R—%. 

^ GT 



(28) 



Notice that when a; — > oo, Z ^ R. 

The dynamic impedance, Z{uj) ~ dV/dl, is easily measured experimentally. It can be determined most readily by 
adding a small A.C. signal to the bias voltage and measuring the transfer function of the detector TF{uj). This is the 
ratio of amplitudes and relative phase between changes in the detector voltage and changes in the bias voltage as a 
function of frequency. Most spectrum analyzers have the capability of measuring the complex ratio between two signals 
as a function of frequency, and can do this simultaneously over the frequency range of interest using a band-limited 
white noise source. Signal averaging allows very precise measurement to be made while remaining in the small-signal 
limit. The dynamic impedance is easily derived from the transfer function using the value of the load resistance and 
making appropriate corrections for stray electrical capacitance or inductance in the circuit. For example, in the bias 
circuit of Fig. El without stray capacitance, the impedance is equal to: Z{ijj) — Rl ■ TF{uj)/{l — TF{llj)). 

It is then possible to determine values for many of the important parameters of the detector by fitting the real and 
imaginary parts of the transfer function by adjusting the thermal and electrical parameters in the expressions given 
in this paper. This is very valuable for diagnosing performance problems or improving the design of detectors. Note 
that when the thermometer temperature coefficient a is positive, the impedance can become infinite. It is then more 
convenient to work with the inverse quantity 1/Z{uj) = dl/dV. 

In the case of a detector whose signal is read out as a voltage change, where the responsivity S is defined as 
S{iLj) — dV (uj) / dW (uj) , we can also write: 



1 (Zo/i?)-l 



21 [Zq/Rl) + 11+ juOTeff 



(29) 



At this point we want to introduce a useful technique for analyzing the response of a bolometer or a microcalorimeter: 
block diagram algebra. This technique is generally used in electrical engineering to analyze feedback systems and it 
is very useful when extending the theory of bolometers and microcalorimeters to more complicated realistic systems. 
The algebra of block diagrams and the language of control theory have been successfully used before in the analysis 
of microcalorimeters and bolometer The basic idea is that a system with transfer function in the frequency 

domain equal to H(u)) is represented by the diagram of Fig.|3| If an input In(uj) is applied to the system the output is 
Out{uj) = H{lo) ■ In{uj). Comphcated systems can always be reduced to the system of Fig. fusing the block diagram 
algebra. Fig. 0] shows some of the common operations that will be used in this paper. The procedure to solve the 
response of a system using the block diagram algebra is then the following: 

• Write the differential equations that define the system response. 

• Convert the equations to the frequency domain and, for each equation define the individual system response 
and the input to that system. 

• Layout the block diagram that describes all the equations together. 

• Use the block diagram algebra to reduce the block diagram to the form of Fig. |21 that represents the system 
response in the frequency domain. 

This representation is particularly useful to deal with feedback systems, i.e. systems where the output is combined 
to the input through a transfer function G{uj) as in Fig. [S| In this case, whenever an external input Inibj) is applied 
to the system the output is: 



Out{uj) 



H{co) 



l + H{uj)G{uj) 



In{Lu) = Hc-i{lu) ■ In{uj), 



(30) 



where Hc-i{uj) is called the closed loop transfer function. 



G l + jfflT 



AT((B) 



X ■ a ■ A,, 



AX(a)) 



FIG. 6; Block diagram representation of a detector. 



Going back to the theory of bolometers and niicrocalorimeters, we can write Ea. ll2l as: 

C^^^^ + G-AT = W~ Getf ■ AT, (31) 
dt 

that, in the frequency domain, becomes: 

juCAT + GAT = W -GetfAT. (32) 

We now want to generate the block diagram describing this equation. The left part represents the response of the 
system that we are analyzing (the output AT as a function of an input power), while the right part represents the 
input to that system. The fact that the input depends on the output AT is a consequence of feedback. The left part 
of the equation represents a low pass system with transfer function 

H{cj) = ^—^ (33) 

G 1 + JUJT 

The input consists of an external input W minus the output itself modified by the transfer function Getf- This is a 
typical feedback system represented by the block diagram of Fig. |H| where we also included the conversion of AT to 
AX. If we now solve the block diagram using the block diagram algebra and Eg. I3UI we obtain: 

A^(t^) = 77— 7; 7— c 7 = -——-. W{uj) (34) 

G + GeTF T Geff l+JUJTeff T 

that is the same expression of Eq. |^ 



II. THE HOT-ELECTRON MODEL 



A first order correction to the standard theory of bolometers and microcalorimeters is the introduction of the 
hot-electron model. The model assumes that the thermal coupling between electrons and lattice in the sensor at low 
temperature is weaker than the coupling between electrons, so that the electric power applied to the electrons rises 
them to a higher temperature than the lattice. This behavior is a known property of metals and has recently been 
quantified in doped silicon so that it affects both TES and semiconductor sensors. The detector can therefore be 
described as composed of two different systems: the electron system and the phonon or lattice system and the two are 
thermally connected by a thermal conductivity Ge-i- We assume for models derived in this paper that the detector 
resistance responds to the temperature of its electron system, and that the Joule power of the bias is dissipated there. 
For economy of presentation, the models derived here assume the input power enters through the absorber phonon 
system, which is then thermally connected to the thermometer lattice, and further to the heat sink through the 
thermal conductivity G (see Fig.O. There are important classes of detectors where signal power is absorbed directly 
in the electron system of the thermometer or absorber, and the primary thermal path to the heat sink could be from 
the absorber lattice or either electron system. In the general case, these all result in different thermal circuits, and 
the block diagrams must be modified accordingly. In the approximation of this section, the phonon system includes 
both the absorber and the phonons in the thermometer (we will discuss later the case of a decoupled absorber). 

In equilibrium with no external power, the electron system is at a higher temperature than the phonon system 
is and the Joule power flows from the electron system to the phonon system and from there to the heat sink. The 
equilibrium temperature of the two systems without any signal power applied can be calculated in a way similar 
to that used for the simple model described in the previous paragraph. As reported in the literature, the thermal 
conductivity between electrons and phonons can be described as a power law of the electron temperature Te [Tol ITU : 



(35) 



P=fR 



'e-l 



w 



'/////////. 

HEAT SINK 

FIG. 7: Thermal sketch of a bolometer or microcalorimeter in the hot electron model 



From the definition of thermal conductivity, we also have: 



dP 
IE 



(36) 



and if we combine the two and integrate from the lattice temperature Ti to the equilibrium electron temperature Tg 
we obtain 



Tf=+i = %±l.p(T,) + rf=+i 

Goe 



(37) 



where we explicitly indicated the dependence of the power P on the electron temperature Tg. The equilibrium 
temperature of the lattice system is still determined by Eq. ^that, in this case, can be written as: 



/3+1 



(3+1 
Go 



P{Te). 



(38) 



Equations 1371 and 1381 represent a system with two variables and T; that can be solved numerically. 

Here we are considering detectors where the external power is absorbed in the phonon system, and the sensitivity 
of the detector can be strongly affected by the reduced sensitivity of to changes in T; introduced by the equilibrium 
difference of these temperatures and non-linear nature of Ge-i- We consider these effects in two steps. A first 
approximation is to assume that the heat capacity of the electron system is negligible. This case can be solved easily, 
and it is sufficient in many cases. We will then derive the general result for Cg ^ 0. 



A. Hot Electron Model with Ce = 



If the electron system heat capacity Ce can be neglected, the dependence of the electron temperature on the 
lattice temperature is simply determined by Eq. |23 When the temperature of the lattice system changes by AT; the 
temperature of the electron system will instantly change by ATg and Eg. 1371 becomes: 



(T, + ATe)^=+i = • P{T, + ATe) + [Ti + 



Goe 



If we subtract Eq. EH from Ea.l^ we obtain 



(Te + AT,f'+' - Tt+^ = ^1±1{P{T, + AT,) - P{T,)) + [Ti + ATif^+^ - t(''+\ 

Goe 



(39) 



(40) 



Assuming that AT, < Tg and AT; < Ti we can expand Eq.^Hlto lowest order in ATe/Te and ATi/Ti, obtaining: 



(1 + + 1)^) . +1 - rf=+i = • AP + (1 + (/3, + 1)^) . Tf^+i - rf=+^ 



Goe 



Ti 



(41) 



1 

G 1 + jtOT 



AT|(a)) 



ATe(Q)) 



AX(a)) 



(a) 



W(m)4 



1 1 



G l + jCOT 



AT|(to) 



^c-IGetf 



X ■ A j_l ■ a ■ A 



AX(ffl) 



(b) 

FIG. 8: Block diagram representation of a detector using the hot electron model with Ce ~ 0. (a) Block diagram as derived 
from the the equations that describe the detector. Notice that the representation of the ETF as acting on the lattice system 
of the sensor is due to the fact that we are assuming that Ce = 0. In general, if Ce 7^ 0, the ETF is an electric effect and acts 
on the electron system, (b) Equivalent representation to highlight the effect of the term A^-i- 



that reduces to: 

Tf= . ATe = — + rf = • AT;. (42) 

We already calculated the change in Joule power AP in Eas.li^ and 1101 that in the hot electron case depends on 
the change in electron temperature ATg: 

AP = -Getf^T,, (43) 

therefore 

Ti" Arj. _ Ge-l{Tl) ^^_Ge-l{Tl) 



ATe = „ ' ^ AT, = , ;^ ' AT, = AT, = Ae_,AT, (44) 

Tt + ^^ Ge-l{Te)+GETF Ge-l{T,)l + ^g^ ^ ^ 



where Ge-i{Te) is the electron-lattice thermal conductivity calculated at the electron temperature and Ge-i{Ti) is 
the electron-lattice thermal conductivity calculated at the lattice temperature, and: 



ATe _ G,_,(T,) J_ 

G 



~ AT " Ge-iiT,) 1 + ■ ^^^^ 



The quantity A^-i is adimensional and represents the temperature sensitivity of the thermometer. When A^-i = 1, 
the thermometer is completely sensitive to temperature changes in the lattice system, when Ae-i = the thermometer 
is completely insensitive to temperature changes in the lattice. 

We now want to represent the detector using the block diagram algebra. The detector behavior is described by 
Eqs.^ll^ and 1441 that in the frequency domain can be written as: 

jioCATi + GATi ^W- GETpATe, (46) 



ATe =Ae_/AT,, 



(47) 



and 



AX ^ ^^ATe 

T, 



(48) 



Converting these three equations in block diagram algebra and connecting the blocks of the algebra together we obtain 
the representation of Fig. |SJi. With some simple algebra, the diagram is equivalent to that of Fig. and considering 



AT,(a)) 



G,_,(T,)Ur^_J 



AT,(Q)) 



FIG. 9: Representation of the electron system in the hot electron model with Ce 7^ 0. 



that Getf oc a, the hot electron model with negligible heat capacity of the electron system is then equivalent to the 
standard model with the substitutions: 

a aeff = Ae-i ■ a (49) 



T^Te (50) 

and therefore the responsivity of the detector becomes 

1 1 Ae-i ■ a - X ■ Atr 

" (G + A,_iGetf) (1 + jWejy) " ^^^^ 

with Tejff — Ci/{G + A^^iGetf), where C; is the lattice heat capacity. 

B. Hot Electron Model with / 

If the heat capacity of the electron system is not negligible, the electron temperature is defined (in analogy to Eq.O 
by: 



Ce^+J^ Ge-l{T')dT' ^P{Te) (52) 

with Ge-i defined by Ea. l35l What we are interested in is the rise of the electron temperature AT^ above equilibrium 
when the lattice temperature rises by AT; . Eq. 122 then becomes 



d{Te + ATe) 



Ge^^^^^—=^+ I Ge-l{r)dT' = P{Te + ATe). (53) 



Subtracting Eq. [S^l from Eq. we obtain 



d(Ar.) , 



Ce^^+y Ge-iiT')dT' - j ^ Ge-iiT')dT' ^AP (54) 

which, using Eg. 1351 and expanding the result to lowest order in ATg/Te and ATi/Ti, becomes 

(j^di^el ^ Q^_^^T^)^T^ ^ G,-i{T,)ATi - Getf AT,, (55) 
dt 

with Getf defined by Eq. 031 The system of Eq. [331 is represented by the block diagram of Fig. |51 with Te_; = 
Ce/Ge-i{Te), which has the solution 

ATeiu;) = A,^i—^. AT,, (56) 

1 +JUJTe 

with Te = Ce/{Ge-i{Te) + G etf) and Ae-i defined by Eq.liH Notice that Eg . reduces to Eq.|33|if Ce = 0. 

The behavior of the lattice system is still regulated by Eqs.ElHandlHl with the substitutions of Gi for C and of Ti 
for T, i.e.: 

Gi^^^ + G-ATi = W + APi. (57) 



W(g))+^ + 

(a) 



1 1 
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FIG. 10: Block diagram representation of a detector using the hot electron model with Ce 7^ 0. (a) Block diagram as derived 
from the equations describing the system. (b),(c) Intermediate steps for the solution of the block diagram representation. 



The power Pi is the power flowing from the electron system to the lattice system through the thermal conductivity 

Pi= I ' Ge-i{T')dT'. (58) 



Therefore 

AP,= / Ge-i{T')dT'- I Ge-i{T')dT', (59) 

Jt^ Jt, 

which, considering the expression of Eg. 1351 for the thermal conductivity and expanding the result to lowest order in 
ATe/Te and ATi/Ti, becomes: 

APi = Ge-i{T,)AT, - Ge-iiTi)ATi. (60) 

Eq. [S7|then becomes: 

a^^^ +G-ATi = W + Ge-i{T,)AT, - G,-i{Ti)ATi. (61) 

Eauations l55l and l61l can be written in the frequency domain as: 

jLoCeAT, + Ge-i{Te)AT, ^ Ge-i{Ti)ATi - GetfAT, (62) 

and 

jijGiATi + G-ATi = W + Ge-i{Te)ATe - Ge-i{Ti)ATi, (63) 



and are represented by the block diagram of Fig. llQb . The diagram can be solved to obtain an analytical expression 
for the detector responsivity. In Figs.llOb and II Oh we show two intermediate steps in the solution of the block diagram 
algebra. The detector responsivity is then equal to: 

Siu)^ ^ ^e-i-X-a-A,^ _ 

GETFA,^i{l+jiUji^) + Gil+ju;n){l+ju;T,) 



1 1 



G l + jtBT 
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ATe(Q)) 



IRa 



Xa-A, 



+ 

j; AX((B) 




IRa 



1 1 


AT(Q>) +^ 


G l + j(BT 





eamp(®) 

+ 




(b) 

FIG. 11: Block diagram representation of the noise in a a detector using the hot electron model(a) and equivalent representation 
for the ideal model (b). Notice that if the output X is a current, the load resistor noise that adds to the output is represented 
by ■ 



with Ti — Ci/G. Notice that in the case of Cg = this expression reduces to Eq. |^ i.e. the hot electron model with 
negligible electron heat capacity, as expected. 

Moreover, in the case of Gqc — > oo, i.e. electrons and phonons can be thermally considered as a single system the 
responsivity becomes: 

y G+Getf J 

This is just the ideal responsivity of a bolometer or microcalorimeter with thermal conductivity G, temperature Tg 
and heat capacity C = Ci + Ce- 

In analogy to Eas. l23lthrough [77l we can also calculate the dynamic impedance of the detector. We can write Eas.l62l 
and 1631 without external power W and explicitly using the symbol AP for the change in Joule power: Converting 
Eos. 1551 and 1611 to the frequency domain we obtain: 

JwCeATe + G,_;(re)ATe = G,_i{Ti)ATi + AP (66) 

and 

jiuCiATi + GATi = Ge-i{T,)AT, - G,-i{Ti)ATi. (67) 
Combining Eqs. □ ESI EH EiandEH we then obtain: 

{G + G,-i{Ti)+ju:Ci)(Ge-.^ . , , 

Z(uj) = R ) U- (68) 



(G + Ge- 


i{Ti) - 
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f jujGe + 
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^JL0Cl)(^Ge- 




\- jioGe - 
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III. NOISE SOURCES 



There are several noise sources that affect the performance of bolometers and microcalorimeters, most of which have 
already been taken into account by Mather in 1982 jlj]. These include the Johnson noise of the sensor, the thermal 
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noise due to the thermal link between the detector and the heat sink (also referred to as phonon noise), the Johnson 
noise of the load resistor used in the bias circuit and the noise of the read-out electronics (amplifier noise). In his 
paper Mather also mentions a 1// noise contribution that seems to be more related to the sensor characteristics. This 
noise has been studied and quantified for silicon implanted thermistors by Han et al. in 1998 T^]. 

The effect of the noise on the detector performance is generally quantified by the Noise Equivalent Power {NEP). 
The NEP corresponds to the power W{u}) that would be necessary as input of the detector to generate an output 
l^X{u)) equal to the output generated by the noise. The NEP is calculated as the ratio between the output AX(a;) 
generated by the noise and the responsivity of the detector S{u)). In the case of bolometers, the NEP directly 
quantifies the limit of the bolometer in detecting a power signal at frequency w. In the case of microcalorimeters the 
NEP is related to the best possible energy resolution of the microcalorimeter by the expression : 

/\Erms = , ^ (69) 



2dui 



NEP2{u) 

Here we want to analyze the effect of the noise on the detector performance in the picture of the hot electron model. 
The introduction of the hot electron model has two main effects: it changes the NEP of the noise sources and it 
introduces a new noise term, that is the thermal noise due to the thermal fluctuations between the lattice and electron 
systems. 



A. Effect of the Hot Electron Model on the Noise 



The different noise contributions affect the detector in different ways. In particular, the thermal noise corresponds 
to a power noise on the lattice system. The Johnson noise is calculated as a voltage fluctuation but can be introduced 
in the model as an electron temperature noise term. The 1// noise is calculated as a fluctuation in the value of the 
resistance but can be described as electron temperature noise term as well. The load resistor noise can be described as a 
noise that adds to the output signal and also generates a Joule power noise on the electron system. The amplifier noise 
adds directly to the output signal. In Fig. Illh . the contribution of the different noise terms in the microcalorimeter 
are shown. The same noise sources in the ideal model scenario are shown in Fig. Illb 9] . Dimensionally, the thermal 
noise Pth is a power spectral density (in units of W/\/Hz), the Johnson noise ej and the load resistor noise brj^ are 
voltage spectral densities (V/y/Hz), the 1/f noise has dimensions Hz'^l"^ and the amplifier noise Camp has 

the dimension of the transducer output X divided by square root of frequency {V/ \J Hz or Aj \J Hz). The thermal 
noise has been calculated quantitatively by Mather in 1982 (assuming diffusive thermal conductivity) and is equal to: 



' ^^"^OTf I '^' -^ I (70) 

JTs kiTi)""^ 

where kb is the Boltzmann constant and k{T') is the function describing the temperature dependence of the thermal 
conductivity of the heat link material. 

The Johnson noise of the sensor resistance is simply described by: 



ej = y/AkbTeR (71) 
The Load resistor noise can be represented by a voltage noise across the detector, equal to: 

R 



CR, = V'^hTsRLj^^-^ (72) 

where we assumed that the electrical circuit is heat sunk at the temperature Tg. This noise adds directly to the 
output signal as a voltage cr,^ or as a current ir^^ = eR^/R, and generates Joule power noise in the electron system 
PR^=2I-eR^ (see Fig. CH). 

The 1/f noise is, by definition, frequency dependent, and it is usually described as a fluctuation in the value of the 
resistance: 
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Solving the block diagram of Fig. II II independently for each noise contribution and using the expression of S(uj) of 
Eq. |M1 we obtain: 

NEPth - Pth (74) 



NEPn, = 1^ + pl^^ (G + G,_KT/) + J^a) (76) 

S[UJ) Ge-l{Tl) 

/A/?\ 7^ 7^/3e / ^0e\ 

NEP,/f = (^— j^^^ (^G(l + ,c.rO(l + ,-c.re-0 + J-G.^J (77) 

iVi^P„„p = (78) 

Notice that the NEP due to the read-out electronics and to the load resistor are the only terms that depend on 
the electro-thermal feedback. Therefore, if these terms are small compared to the other contributions, as it is usually 
the case, the electro-thermal feedback changes the time constant of the detector, but does not affect NEP{lo). 

The expression of the NEP in the case of negligible electron heat capacity is easily derived using Ce = 0. Notice 
that in the limit of Goe oo. Eg s . VP^ and ITSl reduce to the ideal expressions calculated by Mather Kjyj]: 

NEPth = Pth (79) 



NEPj = ^^T^G ( 1 + ju^^ ] . (80) 



Pa ^ \ ^ G 



B. Thermal Noise Due to Hot Electron Decoupling 



The hot electron model also introduces an extra noise term in addition to those just considered. This is due to 
power fluctuations between the lattice and electron system. The magnitude of these fluctuations depends in part on 
the physics of the electron-phonon decoupling. A simple expression appropriate for "radiative" energy transfer was 
calculated by Boyle and Rodger in 1959 



opiin 



Phe = j2hGe^l{Te) ^ ^3 ' . (81) 

A more rigorous expression for electron-phonon decoupling has also been calculated by Golwala et al. in 1997 [T^ . 

Notice that these fluctuations transport power from the lattice system to the electron system and viceversa, therefore 
if a power Phe adds to the electron system, the same power Phe is subtracted from the lattice system. The effect is 
shown in the block diagram of Fig. ^2 Solving the block diagram for the hot electron noise we obtain: 

NEPhe = Phe ^ ^ (l+JLUTi), (82) 

where r; has been previously defined as t/ — Gi/G. This expression does not depend on Ce, therefore is valid also for 
the case Ge = 0. Moreover, if Gqe — oo, this term is zero, as expected. 



w 



G 



'/////////, 

HEAT SINK 

FIG. 12: Thermal sketch of a bolometer or microcalorimeter in the case of absorber decoupling and hot electron model 

IV. ABSORBER DECOUPLING 

Another aspect that may affect the performance of bolometers and microcalorimeters that we want to study is the 
effect of the absorber thermal conductivity. Most of the detectors are built with absorber and sensor as different 
entities connected by epoxy or other material with a thermal conductivity Gq. Depending on the experimental setup, 
there are different configurations that must be used to describe the thermal system. For example, the thermal link to 
the heat sink can be through the absorber or the thermometer and the absorber can be in thermal connection with 
the lattice system (when an electrical insulating material is used) or the electron system (when a conducting material 
is used). 

What we want to analyze here is the case in which the detector is connected to the heat sink through the thermometer 
lattice system and the absorber is connected to the lattice system of the thermometer. In this case, the external power 
hits the absorber, is released to the lattice system and then detected in the electron system (see Fig. lT^ . We assume 
that the absorber has a heat capacity Cq. Notice that the analytical tools that we give here can be easily used to 
quantify the behavior of any other configuration. 

A. Responsivity and Dynamic Impedance 

In equilibrium, with no other power input than the Joule power in the sensor, there is no power flow through the 
thermal link Ga and therefore the temperature of the absorber is equal to the lattice temperature Ta — Ti. If an 
external power W is applied to the absorber, the detector is described in the frequency domain by the set of equations: 



JUjCaATa + GaATa + GaATi 



(83) 



juGiATi + {G + Ga + Ge-i{Ti))ATi = G^AT, + G,.i{T^)AT, 



(84) 



juCeAT, + Ge-i{T,)AT, = Ge^i{Ti)ATi ~ GetfAT, 



(85) 



If we want to build the block diagram associated with these three equations we can consider the left side of the 
equations as the response function of the three systems (absorber, lattice, electrons), and the right side as the input 
to each system. Connecting the three systems gives the block diagram of Fig. El The diagram can be solved to 
obtain the detector responsivity: 



S{u) 



A^^iXaAtr 



iwCa(l+jwr,) + (l+iwre)(l+ia)r„)(G+Ge_,('r!)+i"C',)-G,_i(7'e)A,_,(l+jc^r„) 



(86) 



with Ta = Ca/Ga- Noticc that if Ga oo this expressions reduces to the one without absorber decoupling for a 
detector with lattice heat capacity Gi = Ci -\- Ga- 
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FIG. 13: Block diagram representation of a detector with a finite thermal conductivity between absorber and lattice system. 
We have used the notation Gi = G + Ga + Ge-i{Ti). Notice that this implicitly integrates the heat relief for the lattice system 
provided by the electron and absorber decouling into the lattice response function. This is different from what was done before 
in Fig. 1101 where the heat relief was explicitly reported in the block diagram as a feedback effect. The two descriptions are 
equivalent. We used the implicit description here to compact the block diagram algebra. 



Using Eqs.m53|5Sl and 1831 through we can also calculate the detector dynamic impedance: 

. ^ _ [(G+G,-KrO+3c^CO(l+Jt^r„)+j^C„](G,-KJ^e)+3c^Ce + |f )-G,-,(7^e)-Ge-i(T,)-(l+jt^r„) 
^ [(G+G,_KTi)+jL^CO(l+Jtjr„)+ja;Ca](G,_K'r.)+j<^C.--ef )-Ge-K're)-G,_K'ri)-(l+ja;ra) 



(87) 



B. Noise contribution 



As in the hot electron model of the thermometer, there are two effects introduced by the thermal link between the 
absorber and the lattice system. The first effect is that the response of the detector is different, therefore the NEPs 
due to thermal, Johnson, 1//, load resistor, amplifier, and hot-electron noise are different. The second effect is the 
introduction of an extra noise term due to the power fluctuations between absorber and lattice. 

Fig. 1141 shows the block diagram of the detector with the noise sources evident. As in the hot electron model, the 
noise due to the link between absorber and lattice can be described as a power flow out of the absorber and into 
the lattice or viceversa. This power has same value Pa, but opposite sign at the two ends of the link. Since the 
temperature of absorber and lattice systems are equal, the value of Pa is simply: 



Pa = ^AkbGaTf. 

Solving the block diagram of Fig. 1141 independently for each noise source we obtain: 

NEPth^Pth{l+3L0Ta) 



(89) 



(1 + ]UJTe-l) [jiuCa + G(l + jLUTa)il + JUJTl)) + (1 + jW„) -^^JwC, 



(90) 



1// 

(1 + jujT^-l) {jojCa + G(l + juJTa){l + juTi)) + (1 + j^Ta) ^ jcjCe 



(91) 



NEPr, 



21 e I. 



Siij) Ge-l{Tl) 



[juCa + (1 + JOJTa){G + Ge-l{Tl) + ]uCl)] 



NEPamp = eamp/S{uj) 



(92) 
(93) 



NEPhe = 



he 



Ge-liTi) 



[jujCa + G(l + jujTa){l + jojTi)] 



(94) 



NEPa = PajuJTa 



(95) 



Notice again that if Ga oo these expressions are equal to the hot-electron expressions for a detector with lattice 
heat capacity C/ -I- Ca and the absorber NEP is equal to zero. 



FIG. 14: Block diagram representation ol noise in a detector with a finite thermal conductivity between absorber and lattice 
system. 



V. NON-OHMIC BEHAVIOR OF THE THERMOMETER 



Another effect that may change the performance of a detector is the non-ohmic behavior of the thermometer, i.e., 
the thermometer resistance may not depend only on the thermometer electron temperature, but also on the current 
(or voltage) that is used to readout the temperature change: R = R{Te,I) This effect is particularly strong 

when TES thermometers are used 13]. The responsivity of a detector with non-ohmic thermometer has already 
been calculated by J. C. Mather in 1984 ^3 and its effect on TES microcalorimeters has been studied in detail by 
M. A. Lindeman in 2000 A non-ohmic thermometer can also be easily included in our model. If the resistance 
of the thermometer depends on the readout signal, we can write: 



R R 
dR — —ajdTe + —(3idl 
1 p 1 



(96) 



or equivalently 



R R 
dR = —aydTe + —f^ydV, 



(97) 



where: 



Te dR 
'RdTp 



]_dR 
R~dl 



ay = 



Te dR 
R dTp 



and Pv = 



V dR 
'RdV 



(98) 



Using Eq.^nior Eq.|23is equivalent, and it is always possible to go from one notation to the other using Ohm's laws: 



and av 



ai 



(99) 



The only terms in our model that are affected by the non-ohmic behavior are the electro-thermal feedback term 
Getf and the transducer responsivity Af^. We can calculate them assuming the bias circuit of Fig.|21 



P 



I^R 



AP = 2IRAI 



I^AR, 



(100) 



/ = ^''"^ ^ A/ = Ai?, 



Rl + R 



Rl + R 



(101) 



and 



V = Vuas - IRl ^ av ^ -RlAI. 



(102) 



Using Eq. I^we obtain: 



(103) 
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— ="^i?,+i?(l+/?,)^' 

and 

"T = -"^i?, + i?(i + /?,)—■ 

The model describing a non-ohmic thermometer is therefore identical to that describing a linear one, with the 
substitution: 

ai (106) 



and 



for voltage readout, or 



for current readout. With this substitution in the equations that we derived previously in the paper, it is possible to 
predict both responsivity and noise in the detector. 

We can also use Eq.|^to calculate the dynamic impedance of the detector. In the case of absorber and hot-electron 
decoupling we obtain: 

Z{lo) = i?(/3/ + l) 

[{G+G,^,{T,)+ju>Ci)il+ju>T^)+3ujC^][G,^i{T,)+jL^C, + ,j,^^^j^ (HO) 



This reduces to: 



^ ^ (C+a.-,(T,)...G,)(G (T.).,.C..^) -G._,(T.).G.„,(T,) 



for hot-electron decoupling only, and to: 



1+ Grffe+i) 



Z{lo) = R{Pi + 1) ^ '''Z: ' . (112) 

1 - "St 

for the ideal model. 

We do not know of a rigorous general method for deriving the Johnson noise in a non-ohmic resistor. Nor does 
there seem to be a single definite scheme for determining the net response of the detector to this fundamental thermal 
noise, since it is an internal noise generated in the non-ohmic resistor, and it is not clear how it should itself affect the 
non-ohmicity of the resistor. We are investigating this further, but for the present have assumed that the Johnson 
noise can be represented as a random voltage source with power spectral density AkbTeR in series with the non- 
ohmic resistance, and that the Johnson fluctuations in the source cause the resistance to fluctuate due to the current 
dependence of the resistor. This results in the same suppression of the Johnson noise due to the current dependence 
of the resistance as occurs for external signals and noise if the non-ohmic resistance is expressed as R{Te,I). This 
uncertainty (or dependence on the details of the physics) applies only to the Johnson noise of the sensor. Small-signal 
responsivities to all external sources of signal and noise are unambiguous, so it is only the detector Johnson noise 
contribution to the NEP that is uncertain. 
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FIG. 15: Comparison between the noise from a 6x6 XRS array pixel (courtesy of Caroline K. Stable) and our model. Tbe 
model includes tbe effect of tbe decoupling between hot electrons and phonons in the sensor and between absorber and sensor. 
The noise sources that are included are Johnson noise of the sensor, thermal noise due to the link between detector and heat 
sink, thermal noise due to the link between phonons and electrons in the sensor, thermal noise due to the link between absorber 
and sensor, Johnson noise of the load resistor, 1 // noise, and noise of the readout electronics. 



To verify our results we simulated the performance of an existing microcalorimeter and compared the results with 
data from the detector. We considered a microcalorimeter used in the development phase of the X-Ray Spectrometer 
(XRS) for the Astro-E satellite J^]. The detector that we used for the comparison was part of a 6x6 test array of 
microcalorimeters with silicon implanted thermistors and HgTe absorbers. We chose this detector because the array 
has been studied in great detail and the characteristics of the pixels are well known. 

We first used Eq.|3Zland|3Hlto calculate the expected equilibrium temperature of the detector. We then used Eqs.lHHl 
through EHl to calculated the expected noise spectra. The sum of these can be compared with the measured noise 
spectrum, as shown in Fig. 1151 In the model all the input parameters are fixed to the values measured experimentally. 
The only value that was not available and that has been adjusted during the calculation of the theoretical noise is 
the stray capacitance between gate and source of the FET electronics. The value of 5 pF obtained for the stray 
capacitance is in good agreement with typical values for the FET amplifiers used in the measurement. The agreement 
between the model and the measurement is very good. 

The data set has been acquired at a heat sink temperature of 65 mK. The model predicts an equilibrium temperature 
of 77 mK and, through Eg. 1691 an energy resolution of 8.4 eV, to be compared with the measured values of 78 mK 
and 8.65 eV. The agreement is well within the accuracy of the input parameters in the model and demonstrates the 
power of the model in predicting detector performance. 

Our analytical model has also been compared with a model that uses matrix notation to numerically solve the 
linearized differential equations of the microcalorimeter . The numerical model was developed independently at 
the NASA/Goddard Space Flight Center to predict the performance of more complex detectors Using the same 
parameter values, the agreement between the two is within the numerical error in the implementation of the models 



We have developed an analytical model that predicts the behavior of microcalorimeters and bolometers. The model 
includes the effect of hot electrons in the detector sensor, a thermal decoupling between absorber and sensor, and 
the effect of a non-ohmic thermometer. The model analytically predicts the detector responsivity and expected noise 
under these conditions. The noise sources that are included in the model are the Johnson noise of the sensor, the 
thermal noise due to the link between detector and heat sink, the thermal noise due to the link between phonons 
and electrons in the sensor, the thermal noise due to the link between absorber and sensor, the Johnson noise of the 
load resistor, the 1// noise as thermal noise, and the noise of the readout electronics. A comparison between the 
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predictions of our model and data from a detector developed for the XRS instrument shows good agreement. 

We also described a different way to analyze the performance of bolometers and microcalorimeters, using block 
diagram algebra. The formalism that we introduced can be applied to the description of different detector configura- 
tions. 
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